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\Q . Abstract 

The analysis presented in the recent publication of the CRESST-II results [l[ 
finds a statistically significant excess of registered events over known background 
contributions in the acceptance region and attributes the excess to a possible 
Dark Matter signal, caused by scattering of rather light WIMPs. We propose 
r*] ' a mechanism which explains the excess events with ion sputtering caused by 

<*S || 206 Pb recoils and a-particles from 210 Po decay, combined with realistic surface 

J^ ■ roughness effects. 
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£> | 1. Introduction 

vo' 

[ — ■ The experimental search for Dark Matter has recently entered a new excit- 

l/"~) , ing phase, reaching sufficient sensitivity to directly challenge some of the pro- 

posed theoretical extensions of the Standard Model. One of the Dark Matter 
ff) \ candidates, which has attracted particular attention is the so-called Weakly In- 

^D ■ teracting Massive Particle (WIMP). More than a dozen experiments worldwide 

are currently aiming at direct detection of this class of particles. Current best 
limits on WIMP scattering cross-section have been reported by XENON100 [2[ 
and CDMS-II [3[. No unambiguous evidence for the existence of WIMPs has 
k>( ■ been presented so far,however several experiments (DAMA J4J], for detailed in- 

j_j \ terpretation see also @, more recently CoGeNT @ and CRESST [l() reported 

statistically significant excesses of the number of detected events over expected 
contribution from all known backgrounds, which, under certain assumptions 
could be explained with the WIMP hypothesis. 

Thorough understanding of all backgrounds is clearly necessary for the cor- 
rect interpretation of results - especially when the experiment is not background 
free (and all of the currently running Dark Matter experiments are background 
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limited). Typical sources of backgrounds include: 1) e/7 events 2) neutrons 3) 
a-particles and 4) nuclear recoils, with the last two often caused by a-activity in- 
trinsic to the detector or coming from radon gas and its decay daughters. 222 Rn 
with a half-life of 3.8 days is a particularly problematic background, since it can 
leave its decay products on surfaces, implant them just under the surface or even 
diffuse into materials before it decays. The most problematic radionuclide from 
the 222 Rn chain is 210 Pb with a half-life of 22.3 years and its decay daughter 
210 Po, which itself is an a-emitter with a half-life of 138 days. 

In this article, we focus on the recent result from the CRESST-II Dark Mat- 
ter Search and show that one of the four main backgrounds in the experiment, 
206 Pb recoils from 210 Po decay, was possibly underestimated. Qualitatively, 
taking the 210 Po decay into account with realistic models for surface rough- 
ness (that have been motivated by our work in support of the DEAP [TJ, |§J, |9| 
program), explains the excess events. 

2. CRESST-II detector 

The CRESST-II project (Cryogenic Rare Event Search with Superconduct- 
ing Thermometers) [l|, [l0|, [ll| , located at LNGS in Italy, uses cryogenic detector 
modules which consist of scintillating CaWC>4 crystals operated as bolometers 
(sensitive to phonons) and a nearby but separate cryogenic light detector (a 
silicon wafer with a tungsten thermometer). Simultaneous detection in both 
the "phonon channel" and the "light channel" allows for efficient discrimina- 
tion between electronic, a-particle and nuclear recoils, because of substantial 
differences in light yield. 

In addition to this, the crystal together with the light detector are completely 
surrounded by scintillating and reflecting foil. This is critical for discrimination 
of 210 Po decays, which occur on the inner surface of the reflective housing or on 
the crystal surface. When the recoiling 206 Pb enters the crystal, potentially gen- 
crating a WIMP-likc background event, the escaping 5.3 McV a-particle reaches 
the foil and generates enough light to veto the event. 

As discussed in the CRESST paper, the only non-scintillating surfaces inside 
the detector modules are the small silver-coated bronze clamps holding the 
target crystals, which have been identified as the main source of a's and Pb 
nuclear recoils. This detail is critical for the rest of this discussion. 

3. Analysis from the CRESST-II paper 

The CRESST-II paper extensively covers details of the maximum likelihood 
analysis applied to the final energy spectrum from the entire dataset. The final 
maximum likelihood fit is based on a detailed model including parametrized 
contributions from all backgrounds (as listed earlier) and a possible WIMP 
signal. For each background, a reference region (i.e. range of energies and light 
yields) is defined. A data-driven model or approximated analytical function is 
developed to describe the events in the reference region and eventually used 



to extrapolate the number of background events "leaking" into the acceptance 
region. 

For the e/7-band a Gaussian model is employed, which parametrizes the 
observed energy-dependent light distribution. Because of the high statistics 
available for this class of events the extrapolation to the acceptance region is 
fairly accurate. 

Simple linear functions with parameters constrained by the data are used 
to extrapolate the a-particle spectra, separately for each detector module. The 
spectra are approximately flat over both the reference and acceptance regions 
and the extrapolation is fairly independent on the choice of the reference region. 

A more sophisticated model, which includes hit multiplicities in more than 
one detector, is employed to account for neutrons of cosmogenic origin and for 
neutrons coming from natural radioactivity in the surroundings of the detector. 

Finally, the 206 Pb recoil background spectrum is modelled in a manner simi- 
lar to the case of a-particlcs. A clear peak from the full 206 Pb deposited energy 
of 103 keV is present in the spectrum, with an exponential tail extending to- 
wards lower energies. In the specified reference region between 40 and 90 keV, 
the tail is modelled by an exponential on top of a constant contribution: 
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where Ap h may be module-dependent (to describe the different recoil back- 
ground rates in the individual detector modules). The same spectral decay 
length JSjecay and cons tant background term Cpb are used for all modules. As 
noted by the authors: "the latter quantities are characteristic of the implanta- 
tion profile of a-emitters in the clamps and can thus be assumed to be universal 
if the underlying implantation mechanism is the same for all clamps" . 



The authors used the SRIM package [12( to simulate the energy spectra 
expected for different depth distributions of the 210 Po parent nucleus. The 
simulated spectrum which best fitted the experimental spectrum in the reference 
region corresponded to the depth profile of the 210 Po parent nuclei peaking at 
the clamp surface with a 3nm exponential decay length. It was concluded that 
none of the simulated spectra rose significantly towards lower energies within 
the range of the acceptance region, which then justified a simple exponential 
extrapolation of the Pb background. 

In the next section we will show that the model with a simple exponential tail 
does not correctly describe the 210 Po induced spectrum and leads to a significant 
underestimation of the background contribution to the acceptance region. 



4. Lead recoil background 

4-1- CRESST-II background model 

We start by repeating the key SRIM simulation from Fig. 3 in the CRESST 
paper, which is the main justification for the extrapolation. Decays of 210 Po, 



i.e. a 5.3 MeV a and a 103 keV 206 Pb recoil nucleus, emitted back-to-back, are 
generated in the clamp material according to an exponential depth profile with 
3nm decay length peaking at the surface. 

Silver is the relevant base material for the simulation of 206 Pb recoils com- 
ing from the clamps, as the thickness of the coating (several micrometers) is 
significantly larger than the range of recoils in the material. In addition to 
this, simulations were also performed for bronze. For bronze we simply used a 
"typical" composition provided in the database supplied with SPJM: 89% Cu, 
9% Zn and 2% Pb, with 8.82 g/cm 3 density. For silver, natural abundances and 
10.5 g/cm 3 density were used. 

Since no evidence for material dependence was seen in any of the cases 
simulated for this work and the energy spectra, as well as the fit results, were 
within errors consistent for both silver and bronze, for the sake of clarity we 
only present detailed results for silver. 

Figure [1] compares the published SRIM spectrum (digitized from [l|, Fig. 10]) 
and our results. A similar fit function (Eq. ([I])) is used to describe the spectrum 
in the 40 - 90 keV range, with the resulting E^ = 15.1 ± 1.3 keV for silver 
(or 11.8 ± 1.0 keV for bronze). There is a good agreement between simulated 
curves and the decay constants are consistent with the result from the CRESST 
paper, 13.6 keV, if the uncertainty of about IkeV on the CRESST constant is 
taken into account. 

4-2. Contribution from ion sputtering 

A lot of insight can be gained by considering individual contributions to the 
spectrum in Fig. [T] Apart from the obvious contribution from the polonium 
decay products ( 206 Pb ion and a-particle), there is also a significant fraction of 
energy carried into the detector by sputtering products. 

In the sputtering process atoms are ejected from a base material due to 
bombardment of the target by energetic particles. Sputtering is driven by mo- 
mentum exchange between the ions and atoms in the material. The incident ions 
set off collision cascades in the base material. When such cascades recoil and 
reach the surface with an energy exceeding the surface binding energy (typically 
of the order of a few eV), an atom can be ejected. 

An approximated analytical model for describing sputtering in the cascade 
regime for amorphous flat base materials has been developed by Thompson |13l | 
and a more sophisticated algorithm is included in SRIM. The maximum recoil 
energy available to sputtered atoms is given as AE±, where A = AM\M<2,I{M\ + 
M2) 2 , Mi and M^ are the masses of the incoming and sputtered ions, and E\ 
is the energy of the incoming ion. For silver this gives 0.9-Ei, i.e. up to 93kcV, 
for copper, the dominant component of bronze, up to to 75 keV. 

Depending on the kinematics, the distance of the collision point from the 
surface, and the stopping power, a significant fraction of the base material atoms 
can be kicked off the material with a total kinetic energy, typically distributed 
among several ions, exceeding lOkeV. The model predicts 1/E 2 energy depen- 
dence for the spectrum of sputtered atoms, in the region between the surface 
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Figure 1: Energy spectrum of 206 Pb recoils calculated with the SRIM package 
for a depth distribution peaking at the surface with a 3nm exponential decay 
length: (shaded black histogram) plot digitized from Fig. 10 in the CRESST 
paper, (dark blue points) this work, with silver as the clamp material. The red 
curve is a fit of Eq. (fTJ) to the simulated spectrum in the energy range of the 
reference region from 40 to 90kcV (solid line), extrapolated to the acceptance 
region (dashed line). 
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Figure 2: Main components to the background induced by the clamps: (black 
solid line) contributions from 206 Pb recoils alone, (red dashed line) contributions 
from both 206 Pb and sputtered ions, summed up for each decay, and (blue 
points) the full spectrum, i.e. including also the contribution from a-particlcs 
(equivalent to Fig. [T]). 



binding energy, E^ (typically of the order of a few eV) and the maximum recoil 
energy h.E\ . 

Figure [5] shows the energy spectrum broken down into its main components. 
It is evident that taking into account the energy of sputtered atoms introduces 
a slope increasing towards low energies. A significant energy share carried by 
sputtering products at the low-energy end of the spectrum is an effect that the 
CRESST collaboration has investigated both experimentally and by means of 
SRIM simulations 
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Interestingly, the enhancement at low energies is not seen in the simulated 
full spectrum, as the low-energy events dominated by the ion sputtering prod- 
ucts are typically accompanied by a high-energy a-particle entering the detector, 
as will be discussed later. This effectively moves the events out of the region 
of interest towards higher energies (and also towards higher light yields, as the 
recoil scintillation is significantly more quenched than the alpha scintillation). 

4-3. Surface roughness 

As concluded in the previous section, the flat shape of the spectrum of 210 Po 
decay products at low energies is an outcome of a specific balance between 



sputtered atoms and a's entering the detector. So far, it has been assumed that 
the clamp surface is perfectly smooth. Since ranges of both types of particles 
in silver vary dramatically (~ 10 microns for a's vs. 16 nm for recoils) it is easy 
to imagine that taking realistic surface roughness into account will affect them 
differently and break that balance, with some consequences for the acceptance 
region. 

SRIM does not allow for treatment of arbitrary geometries and rough sur- 
faces, contrary to Geant4 [lfj, a package for simulation of the passage of parti- 
cles through matter, developed by the high-energy physics community. Geant4 
is a flexible framework, with increasing usefulness also for low-energy physics 
applications. In particular, one of its extensions, available as one of its extended 
examples (extended/electromagnetic/TestEm7) from within the standard distri- 
bution, contains all physics relevant for multiple interatomic and alpha scatter- 
ing in lOkeV - lOMeV energy range [17|, |lj|. The code is actually meant "to 
bridge the gap between the effective handling of low-energy processes in sim- 
ple geometries provided by SRIM and the very general framework for nuclear 
events and complex geometries already available in Geant4" (see Ref. [18(). It 
has been extensively benchmarked against SRIM with respect to nuclear strag- 
gling and implantation problems, as well as backscattering, with satisfactory 
agreement. 

A simplified geometry of the problematic clamp region was implemented as 
a bronze or silver volume with a 10 x 10 micron rectangular area and a few 
micrometers thick, located just next to a "detector" volume of a larger area. 
Decays of 210 Po are generated in the material (of the same composition as 
in SRIM) according to the same exponential depth distribution, and then all 
decay products and secondary cascades are tracked and registered when at the 
boundary of the detector volume. 

Initially the simulation was repeated for a flat clamp surface, with results in 
good agreement to SRIM, shown in Fig. [3l 

As for the actual surface roughness, a simplified model was used, which 
assumed perfect sine-like waviness in one direction. Directionality of surface 
roughness is, in fact, a frequent outcome of machining or polishing. The surface 
geometry was implemented using a standard Geant4 volume (G4Polyhedra), 
which forced us to approximate the sine-like waviness with finite length linear 
sections, as shown in Fig. He] 

Generally, a completely arbitrary surface profile could be implemented in 
a similar way. It would also be possible to use a 2D surface profile, directly 
measured e.g. with a profiler or an atomic force microscope (AFM). 

Simulations were run for some arbitrarily chosen values parametrizing the 
roughness, shown as cases #1 - #4 in Table [TJ where "waviness", T, is simply 
the period of the sine function describing the surface and the average roughness, 
R„, is defined as: 



R a = j, 



. 2tt 
A ■ sin — X 
T 



dx = — . 



More realistic sets of values are studied in cases #5 - #8. Based on liter- 
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Figure 3: Geant4 simulation results for a perfectly smooth silver surface: (blue 
points) full energy spectrum from all contributions, (red dashed histogram) con- 
tributions only from both 206 Pb and sputtered ions, summed up for each decay, 
compared to SRIM simulation result digitized from Fig. 10 in the CRESST pa- 
per (shaded black histogram) . The red curve is a fit of Eq. ([T]) to the simulated 
blue spectrum in the energy range of the reference region from 40 to 90keV 
(solid line), extrapolated to the acceptance region (dashed line). 
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Figure 4: Zoomed-in cross section of Geant4 geometries showing relevant event 
types: (a) a typical sputtering-dominated event on a flat surface; Cu recoils 
carrying 10 - 40 keV of energy, with a contribution from a high-energy a (which 
will eventually reach the detector and shift the event out of the acceptance 
region) , (b) a cascade of sputtered Cu ions reaching the detector (total energy 
above 40 keV), resulting in an event with no contribution from a and (c) an 
example of a rough surface (case #7), with a part of the cascade stopped by 
surface irregularities, which effectively shifts the event towards lower energies, 
a's are shown in red, 206 Pb in blue and copper recoil nuclei in green. Electrons 
are not shown. Particles are "detected" and killed after reaching the surface on 
the left side of the plot (which extends much beyond the plot boundaries). 
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Figure 5: Geant4 simulated spectra for all respective "rough" cases from Ta- 
ble [1] Plots contain: (blue points) full energy spectrum from all contributions, 
(black solid line) energy from 206 Pb recoils only (no sputtering and no alphas 
included) . The red curve is a fit of Eq. ((T|) to the simulated blue spectrum in the 
energy range of the reference region from 40 to 90 keV (solid line) , extrapolated 
to the acceptance region (dashed line). 



ature data found on the surface roughness of an "apparently smooth" bronze 
surface [19( and the roughness of electroplated silver coatings [20J, [21| , we ran 
simulations for the average roughness of 0.07 micrometers (2.7 /zinch), which 
was reported for bronze and is also consistent with the smoothest silver coat- 
ing reported in Ref. [20J, Fig. 5]. The relevant range of the waviness parameter 
(0.1-5 micrometers) was specified based on AFM surface scans from 20ll21[. 

As mentioned earlier, there is no statistically significant difference in results 
between bronze and silver. Qualitatively, all but one simulated configuration 
show the same effect, namely, the excess of events at lower energies, when com- 
pared to the "flat" surface model. Resulting spectra for all rough cases from Ta- 
ble Q] are shown in Fig. [SJ 

The only exception from the rule is case #5, where the low energy contribu- 
tion disappears, making the result similar to that of a perfectly smooth surface, 
i.e. case #0, discussed earlier. Some non-trivial dependence on the A/T ratio 
is not a surprise, since smoother surfaces will have smaller A/T ratios. This 
gives an idea on what kind of a surface finish would be required to mitigate 
backgrounds of this type. 



10 







Surface Roughness 


Silver 


# 


A [fJ-m] 


T [fim] 


A/T 


R a \p.m] 


decay L J 


„Pb/ Pb 















12.6 ±1.0 


2.46 


1 

2 
3 

4 


1 
1 

10 
0.020 


5.0 

2.5 

5.0 

0.005 


0.2 
0.4 

2 

4 


0.64 

0.64 

6.4 

0.013 


7.8 ±0.9 
6.2 ±0.9 
4.7 ±1.3 
9.5 ±1.4 


1.29 
1.02 
0.89 
1.29 


5 
6 
7 
8 


0.11 
0.11 
0.11 
0.11 


5.0 
1.0 
0.5 
0.1 


0.022 
0.11 
0.22 
1.1 


0.07 
0.07 
0.07 
0.07 


11.4 ±0.9 
7.8 ±0.7 
9.1 ±1.0 
6.5 ±1.4 


2.45 
1.67 
1.31 
1.1 



Table 1: Roughness for all cases used in Geant4 simulations: amplitude (A), 
waviness (T), their ratio, and average surface roughness (R a )- Two last columns 
summarize the simulation results and contain, respectively: exponential decay 
constant fitted to the reference region (E^ ), and the ratio of number of 
counts in the reference region to the number of counts in the acceptance region 
(nftf/n^.) for silver. Statistically consistent results obtained for bronze are 
omitted. Case #0 corresponds to data shown in Fig. [3] 



We did not sec any similar effect of difference between a smooth and a 
rough surface for a's generated deeper in the clamps, i.e. uniformly from a 
15 micrometer thick surface layer. This can be explained by the range of a's in 
bronze or silver, which is much larger than the range of 206 Pb recoils and also 
larger than the typical roughness scale. 



4-4- Light- quenching factors 

So far we have not considered the consequences of light-quenching factor 
(QF) difference between 206 Pb and the sputtered ions. The QF is defined as 
a light output of a particle, relative to the light output of a 7 of the same 
deposited energy. Nuclear recoil quenching factors in CaWC-4 decrease with 
increasing atomic mass of an ion 22|, [23J ■ CRESST [1| reported values for Ca 
and W of 6.4% and 3.9%, respectively, and 1.4% for Pb (all values with about 
10% relative uncertainty). Based on 23J, Fig. 1] one can infer the QF value for 
silver of around 3%. 

Therefore, the low energy part of the spectrum, dominated by the sputtered 
silver ions, should have about 2-3 times larger light yield than Pb ions. This 
is roughly consistent with Figure 12 in the CRESST paper [l(, where the light 
yield distribution of a possible WIMP signal, corresponding to the likelihood 
maximum Ml, is centered around 5%. 



5. Full model and the CRESST-II results 

Below, in Fig. [6j we present a direct comparison of the low energy spectrum 
measured by CRESST with the spectra obtained from Geant4 simulations 
taking the surface roughness effects into account. 
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Figure 6: Comparison of CRESST data (red points) with CRESST's clamp 
background model (black curve) and our model including surface roughness 
(blue shaded histogram). The CRESST data was digitized from Figures 9 and 
11 in Ref. [l|. Our Geant4 simulation of the clamp background spectrum 
corresponds to case #7 (see text) and was normalized to give the same integrated 
number of 89 counts in the reference region as the CRESST data. The CRESST 
background model is based on a log-likelihood fit of Eq. (fT]) to the data in 
the reference region (40 to 90keV, shown as solid line). Extrapolation of the 
fitted function to the acceptance region (10 - 40keV, shown as dashed line), is 
consistent with what was used in the CRESST analysis to estimate the clamp 
background contribution. 
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The exponential fit to the closest matching spectrum, which corresponds 
to case #7, yields 9.1 ± 1.0 keV for the exponential decay constant. As the 
uncertainty on the slope parameter of the exponential fit to the reference region 
has not been included in the CRESST paper, we repeated the fit, confirming 
the published values, in particular, E^ = 13.72 ± 6.19 keV. It is, thus, within 
error consistent with the slope of the "full cascade" spectrum. 

We note, that the full energy 206 Pb recoil peak at 103 keV from the simu- 
lation has lower intensity, when compared with the CRESST spectrum. Out 
of many possible explanations, perhaps the most likely ones are the surface 
component of 210 Po contamination exceeding the assumed exponential depth 
distribution or surface contamination of the CaW04 crystal directly under the 
clamps. 

However, most importantly, we find a fair agreement between the integrated 
number of counts in the acceptance region between our model (68 counts or 
53 counts, based on case #7 or case #6, respectively) and the CRESST data, 
i.e. 67 counts from all backgrounds (which includes very approximately 40 - 45 
counts attributed in the global likelihood fit to the WIMP signal and the 206 Pb 
recoils taken together). 

Based on our qualitative study, the clamp background is significantly larger 
than estimated in the CRESST background model, and might be the domi- 
nant contribution to the low-energy part of the spectrum. It remains unclear, 
however, whether it is sufficient to entirely explain the excess of events, which 
necessitates further, better informed, analysis. 

In addition to this, we note that it is likely that the same class of low- 
energy events correlated with 206 Pb contamination of the clamps, contributed 
significantly to the earlier spectra published by CRESST as Figures 3 and 8 in 



Ref. }24j . which exhibit similar event excesses at lower energies (attributed by 
the authors entirely to neutrons, for the former, and e/7-band leakage, for the 
latter). 

6. Summary 

We have found a significant contribution to the number of events in the ac- 
ceptance region in CRESST-II data coming from cascades of secondary nuclear 
recoils caused by recoiling 103 keV 206 Pb nucleus from decays of 210 Po present 
on and under the surface of clamps, used to hold the target. The size of this con- 
tribution becomes significant only when a realistic model of surface roughness 
is included in the calculation. 

This is a new type of background which to our knowledge has not been 
previously considered in the context of ultra- low background experiments. It is, 
however, fairly easy to model using some commonly used simulation tools, such 
as Geant4. The study itself has been motivated by similar surface roughness 
effects discovered in data from DEAP-1 |[, a prototype for the DEAP-3600 @] 
dark matter detector with a single-phase liquid argon target. 

Based on the published data, we have qualitatively estimated the expected 
number of this class of events in the CRESST-II dataset, and found it ap- 
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proximately consistent with the excess of events in the acceptance region. The 
expected light yield of sputtered silver ions is also consistent with the light yield 
distribution of these events. Although a more sophisticated study including a 
corrected maximum likelihood analysis is necessary to reach the final conclusion, 
and also specific input regarding the clamp composition, nature of contamina- 
tion and surface roughness would be useful, we interpret this as an indication 
that the 206 Pb recoils could be responsible for the dominant part of the signal 
observed by CRESST-ff . 
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